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Summary 


The process of posterior evaluation is often the most demanding part of data analysis 
pipelines, both in human and CPU hours. We propose a new approach to the 
problem. 


@ The likelihood can be factorized by dependence on either intrinsic or 
extrinsic parameters. 

@ Using information from the detection-pipeline trigger, intrinsic and extrinsic 
samples are drawn separately and at once. 

© The respective components (e.g. waveforms, antennae response) are then 
computed. 

@ The likelihood is evaluated at all combinations of intrinsic and extrinsic 
parameters, by using algebraic operations (i.e. "dot” product), as in the warp 
and weft of a fabric. 


This has implications for parameter estimation; finding best-fit waveforms (con- 
sistency tests); and allowing an optimal search test statistics that includes the 
effects of precession and higher emission modes. 


The Factorization 


e We separate the 15 parameters of binary mergers [1] to intrinsic and extrinsic: 
eine a M1, Mo, $1, 52, L, Gere a a, 0, w, is Q, D 


e The waveform can then be decomposed by modes m, polarizations p, and functions 
of either intrinsic or extrinsic parameters: 
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e where F, is the antenna response to the p polarization, and Do, to and gp are some 
arbitrary values. 


e Employ relative binning [2]: we can pre-compute weights wen Wes (using 
a single reference 6;,;.), then use frequencies f, of length O(10°) instead of f of 


length O(10°). 
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Sampling Procedure 


e Samples drawn using importance / rejection sampling, given the detection trigger. 
Employ quasi Monte Carlo for better convergence. 


e L) and ¢ can be marginalized over. 
e No MCMC / nested sampling required. 


e Current benchmarks: N.,, = 107, N= 10° (see Figure 1). Include 4 modes, 2 
polarization, ~ 400 frequencies for ~ 101° operations. After optimization, we aim 
for a run-time of ~ 1 — 10 seconds. 
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Implication for Search Pipelines 


The optimal test statistic (the evidence ratio [3]) can be computed with samples from 
a single-detector posterior, along with mock-data produced cheaply by timeslides of 
the other detectors, for the background (null hypothesis Ho) and foreground (signal 
hypothesis H,). For example, when three detectors are available: 
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Timeslide detectors (2,3) by FFT for mock data 
(d|h(8)) — (d]h()) + (h(a'™"°") |B) 
Foreground distribution by injections 


of waveforms into mock data 


e The test statistic includes higher modes and precession, and is coherent between 
detectors. 


e From a single event, the probability ratio Pr(Z|H,)/ Pr(Z|Ho) is obtained. 
e The method demands a fast and reliable way to produce samples from a posterior. 


e Comes as a refined search, after a coarser search pipeline. 
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Figure 1: Convergence test of a signal 
injected in Livingston, for fixed 

Neo & 10° and different N. n+ 


Figure 2: Declared event 

GW190623 183652. High 
Pr(Z|H1)/ Pr(Z|Hpo) is in agreement 
with reported IFAR [4, 5]. Sampling 
was performed on Livingston, 
timeslides performed on Hanford. 


Astrophysical Implication 


The optimal test statistic will improve the sensitivity to binary mergers exhibiting 


precession or higher modes. The discovery of precessing binaries would suggest 
dynamical formation, while the lack thereof would hint towards isolated co-evolution 


[6]. 
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